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A note on the method of steepest descents with a remark 
on T. Ljunggren’s paper “Contributions to the theory of 
diffraction of electromagnetic waves by spherical particles’’! 


By PER-OLov LOWDIN 


In his theory of diffraction of electromagnetic waves with the wave-length 2 
around a sphere with the radius a and relative refraction index n, Ljunggren 
evaluates a large number of integrals of the form 


I (0) = fa(r) ede (1) 
b 


by means of the method of steepest descents for large values of the parameter 
0 =2za/j. In this method2, in its pure form, the slowly varying functions. 
G(t) and g(t) are developed as Taylor polynomials around certain decisive 
points, namely the limits of integration and the saddle points ty, where at 
least G’(t)) = 0. Here we wish to point out that, if the expression for a saddle 
point has a complicated character, depending e.g. on one or two parameters, 
it can sometimes be convenient to develop G and g around another fix-point 
T with simpler properties, which is chosen in the immediate neighbourhood of 
the saddle point under consideration.* If this saddle point has the multiplicity 
m, the approximate polynomial for (rt) has still to be of at least the degree 
(m+ 1), and one has further to choose a corresponding integration contour. 
The slight loss of accuracy connected with this modification is often without 
any importance and can, if desired, be compensated by considering even higher 
terms in the asymptotic expansion of the integral (1). 


1 T. Lsunacren, Arkiv f. mat., astr. o. fysik 36A, n:o 14 (1948); Arkiv f. fysik 1, n:o 1 
(1949); in the following quoted as LI and L II, respectively. 

* Developed by B. Rremann, Ges. Math. Werke, p. 409 (1876); Lord Kartvry, Phil. Mag. 
23, p. 252 (1887); P. Dresyz, Math. Ann. 67, p. 535 (1909); Munch. Sitz. Ber. 40, n:o 5 
(1910); L. Brittovurn, Ann. d. l’Ecole Norm. Sup. 33, p. 17 (1916); E. Hoaner, Thesis, 
Uppsala 1925; F. PotzaczeK, Ann. Phys. 2, p. 991 (1929); H. G. BarmrwaLp, Ann. Phys. 6, 
p. 295 (1930); V. p. Corput, Comp. Math. 3, p. 328 (1936); and others. 

° A variation of this idea has previously been employed e.g. in the case of two closely 
situated simple saddle points, where it has been recommended to develop around the point 
Too, Where F’’ (T)9) = 0. The integral can then be expressed in terms of Airy integrals. Cf. 
J. W. NicHotson, Phil. Mag. 19, p. 247 (1910); L. Briztoutrn, loc. cit.; E. Hoaner, loc. 
cit.; F. Porwtaczex, loc. cit.; H. Jerrreys, Phil. Mag. 33, p. 451 (1942); T. LauNGGREN, 
Wate pss on (1948). 
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Let us consider the case when a simple saddle point is situated near to one 
of the limits of integration!, say t = c¢. Developing G and g around the point 
T) = c, we obtain easily the following expression for the main contribution from 
these two decisive points to the integral (1): 


’ 2 RM Radel : 
I (9) = seco V/ mle a ge ewe —2ian(a)}, 
z Cc c 
| oo (2) 


eat ey aoe =o:  f2EUHe day, 
o-- RV se 1 | 
0 


Since the function 7(«) can directly be expressed in the ordinary Fresnel inte- 
grals, the contribution from a simple saddle point near to a limit is essentially 
characterized by these integrals. ? 

The purpose of this note is to show that some of the calculations carried 
out by Launccren (LII, p. 3-8) can be essentially simplified by using the 
formula (2). Let us investigate Ljunggren’s integral 


I(e, n, O) = 


7 


nt 
2) 2 
+ co . ° 
a 4n¢e? (1 : =) S = 1" [av | = Ty SIN Ty, COS Ty {oie GD Hake GD) dt, 
0 0 


1 2 her may (n sin tT, + sin T,)? 


ISD) 


== COS ty" (— 20N — 274 + 27, +. O Cos PY) + 23in 7, — 2% sin 2 


1 
COS T, = — COS 74, 


(L II, formulae 5 and 6) for large values of 9. The main contribution from 
the series in N arises from the term N = 0, where the inner integral has a 
simple saddle point in the neighbourhood of the upper limit of integration 
Tt = 2/2 (above or below this point). Using Brillouin’s method, Ljunggren 
develops his functions G and g around the rather complicated saddle points, 
and the calculations can then be performed only with a considerable amount 
of labour. His final expressions (L II, 25 and 25’) are complicated, and, in 
order to get a survey of his result, Ljunggren then summed the infinite series 
involving J 9 by numerical summation for some typical values of the variable. 
Let us here first consider the inner integral for N = 0: 


x 
2 
sin? Ty SIN T, COS Ty 
Linner 7 ae sna 


(n sin OE ane Ne ee (4) 
‘2 


* Cf. the treatment of the same problem given by BriLLoury, loc. cit. 
2 Cf. BapRwaxp, loc. cit. 
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Developing the functions G and g around the upper limit t, = 2/2 with 
Tt = 2/2, we get : 


ge = 0 isis eee 
pains y (n + 1)2 (5) 
G, =— 2(n—1), Ge = + O cos ¥, qe = 22, | 
and ' 
a=+k cos ¥, (6) 
with the new notation 
: 1 ae 
ee : 7 
. 98 e|/ m= 1 MO 
According to the simplified formula (2), we then directly obtain 
Tinner = 2k e—2%e—1) (2 — Ick cos W-n(k cos YW) + 
2p (7% — ih) (ne 31)? 
+ 2ik cos ¥-n(— k cos ¥)} = 
| ; ue (8) 
s un —2io(n—1))1 _ | i? on (9 Y 
Doig Win ine \ ae | e'“ sin (2k u cos He) ter 
In order to carry out the outer integration, we introduce the function 
n)2 +o 
x --{ h — 2k cos Je fj sin (2 ku cos vy du) av = 
0 0 (9) 


=1—2k [een ekudu =1—(1—¢'”) =er*¥, 
it 


where the integral in the second line has been transformed by changing the 
path of integration and developing the Bessel function J, in a power series. 
Hence we get for (3) the final result 
OF n> 
1 @) = 2 “eto (1 4 ey, 10 
(0, ”, ) vee 3 (n —1)(n + 12% ( ) 
where Ljunggren’s function y now can be expressed by the simple analytical 
formula 


—i- Ss 06 
ore, 2 (11) 


It is remarkable that, for very small angles 0, Ljunggren has obtained the 
approximate formula (L II, 27’) 
n 


EE 
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and that his numerical results according to (L II, 25’) indicate the validity of | 
the relation |y| ~ 1 (L II, 48 and 49). a 

Finally, we will give a comparison between Ljunggren’ s numerical values of 
y and the values obtained from eq. (11). It is suitable to tabulate x as a . 
function of the variable 


n 


ear yrees i, 18) 


(cf. the slightly different choice oF B in LII 21), and the following table gives 
a survey of x for some types values of this variable: 


ie = e—2%B a4 4 
B Analytically Numerically 
according to (11) according to L II 25’. 


0.1 0.980—0.199 2 0.982—0.199 4 

0.3 0.825 —0.565 7 0.826 — 0.567 7 

0.5 0.540—0.841 7 0.540—0.841 2 7 
0.7 0.170—0.985 7 0.169—0.985 7 

0.9 —0.227—0.974 4 —0.231—0.971 @ 

el — 0.589 — 0.808 2 — 0.595 — 0.795 4 

1.3 — 0.857—0.516 7 — 0.864—0.490 7 

1.5 —0.990—0.141 4 — 01907 05007 aca 

1.6 — 0.998 + 0.058 7 —0.997+0.116 2 | 


The agreement is excellent, and each of the results may be considered ag 
a check of the other. Consequently the small loss of accuracy associated with 
the use of the simplified formula (2) must here be without any practical im- 
portance; on the other hand we have gained a considerable simplification of 
the final result. 

This remark was first given by the author on March 19th, 1949, in connec- 
tion with the public discussion of T. Ljunggren’s thesis for the degree Doctor of 
Philosophy. I am indebted to Dr. LiunGcGREN for placing his numerical ma- 
terial at my disposal. 


Institute of Mechanics and Mathematical Physics, University of Uppsala, 
Uppsala. 


Tryckt den 17 oktober 1950 


Uppsala 1950, Almqvist & Wiksells Boktryckeri AB 


370 


